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Abstract 

The hyperscaling property implies that spatially isotropic critical quantum states in d spatial dimensions 
have a specific heat which scales with temperature as , and an optical conductivity which scales with 
frequency as for a; S> T, where z is the dynamic critical exponent. We examine the spin-density- 

wave critical fixed point of metals in d = 2 found by Sur and Lee (Phys. Rev. B 91 , 125136 (2015)) in 
an expansion in e = 3 — d. We find that the contributions of the “hot spots” on the Fermi surface to the 
optical conductivity and specific heat obey hyperscaling (up to logarithms), and agree with the results of 
the large N analysis of the optical conductivity by Hartnoll et al. (Phys. Rev. 84 , 125115 (2011)). With 
a small bare velocity of the boson associated with the spin density wave order, there is an intermediate 
energy regime where hyperscaling is violated with d —>■ dj, where dt = 1 is the number of dimensions 
transverse to the Fermi surface. We also present a Boltzmann equation analysis which indicates that the 
hot spot contribution to the DC conductivity has the same scaling as the optical conductivity, with T 
replacing uj. 
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I. INTRODUCTION 


The anomalous properties of the ‘strange metal’ phase of the cuprates, and other correlated 
electron compounds, have remained a long-standing challenge to quantum many-body theory. 
Strange metals are states of quantum matter whose density can be continuously varied by an 
external chemical potential at zero temperature, but unlike in a Fermi liquid, there are no long- 
lived quasiparticle excitations. It is generally believed that strange metals should be described by 
a strongly-coupled quantum-critical theory [1], but such a proposal immediately faces an obstacle. 
Almost all strongly-coupled quantum-critical states, including all conformal held theories, obey 
the ‘hyperscaling’ property [2]: this implies that the specihc heat, Uy, and the conductivity, a, 
scale as 

; a(a; > T) ~ (1.1) 

where u is frequency, T is temperature, d is the spatial dimension, and z is the dynamic critical 
exponent; we will refer to the conductivity in the w 3> T regime above as the optical conductivity. 
In the important spatial dimension of d = 2, this immediately implies that the optical conductivity 
should be frequency independent, which contradicts the ~ behavior observed in the cuprates 

[3, 4]. 

The scaling arguments can also be naively extended to the DC conductivity, which would then 
imply that a{u T) In d = 2, this contradicts the widely observed ‘linear-in- 

T resistivity’, a ~ T~^. However, DC transport co-efficients are sensitive to constraints from 
momentum conservation, and so the naive application of hyperscaling to DC transport is often 
not valid [5-14]. But this sensitivity does not extend to the optical conductivity, and so the 
observations of Ref. 3 and 4 are the stronger challenge to the hyperscaling property. 

There is a much-studied [15-29] strongly-coupled quantum-critical point which violates hyper- 
scaling: this is the critical point to the onset of Ising-nematic order in a metal in d = 2. A 
closely-related critical theory applies to a d = 2 metal coupled to an Abelian or non-Abelian gauge 
held. We write the properties of the Ising-nematic theory in a suggestive form similar to Eq. (1.1) 

Cy ~ ; a{u > T) ~ with hyperscaling violation, (1.2) 

where z = 3/2 is the ‘fermionic’ dynamic critical exponent (in the notation of Ref. 24). For 
the specihc heat, the hyperscaling-violating dimensionality d^ = 1 has been connected to the 
number of dimensions transverse to the Fermi surface [27, 30]. This value of dt also happens to 
yield the correct behavior of the optical conductivity in Eq. (1.2), although the existing [19, 31] 
physical interpretations of this result are diherent. It is also notable that a ~ is close to the 

experimental observations [3, 4]. 

The above violation of hyperscaling is in a theory with a ‘critical Fermi surface’. On the other 
hand, theories with Dirac fermions, which are gapless only at points in the Brillouin zone do obey 
hyperscaling. 
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FIG. 1. (a) Hot spot geometry, labelling conventions, and choice of x, y-coordinate system in the Brillouin 
zone of the two-dimensional square lattice in which the fermions move. The boundary of the blue area 
denotes the Fermi surface separating the filled particle-like states (blue area) from the the hole-like states 
(white area). Qaf = ("FiVt) is the (commensurate) antiferromagnetic ordering vector that intersects the 
Fermi surface at 4 pairs of hot spots, (b) (Top) Fermi surface patches from a hot spot pair connected 
by Qaf centered at a common origin in momentum space. The two light colored regions are the regions 
occupied by fermions at the two hot spots of the pair respectively, the dark colored region is occupied by 
fermions at both hot spots, and the white region is unoccupied. The arrows perpendicular to the Fermi 
surfaces denote the directions of the Fermi velocities. (Bottom) Under the RG flow, the Fermi surfaces are 
deformed as shown at the strange metal fixed point, and as indicated in Eq. (1.3). The Fermi velocities 
are exactly antiparallel only at the hot spot (k = 0). 


Our interest in this paper is the onset of spin density wave order in two-dimensional metals, 
whose critical theory is described by isolated points called ‘hot spots’ which are connected to a 
gapless Fermi surface (see Fig. 1). This transition is therefore intermediate between the critical 
Fermi surface and critical Fermi point cases. Its field theory [32] has a bosonic order parameter 0 
coupled to fermionic excitations at 4 pairs of hot spots around the Fermi surface. 

In a large N analysis of such a field theory, it was found [32, 33] that at the two loop level that 
the Fermi surfaces near the hot spots became asymptotically nested at low energies. In terms of 
momenta k^, ky measuring deviations from the hot spots, the Fermi surface is given by (see the 
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bottom panel of Fig. lb) 


( 1 , 3 ) 


=*= 111 ( 1 / 14 , 1 )' 


The optical conductivity of the hot spots was computed by Hartnoll et al. [31] in a Eliashberg 
framework, and they found (at variance with an earlier treatment by Abanov et al. [34], and that 
in Ref. 35) a hot spot contribution (t{u) ~ where the exponent rg > 0 was determined by the 
angle between the Fermi surfaces at the hot spots. For the asympotically nested Fermi surface in 
Eq. (1.3), it was found [31] that rg —)■ 0, indicating that the optical conductivity is a constant (up 
to logarithms), and so obeys hyperscaling as in Eq. (1.1) in d = 2. 

This paper will re-examine these issues using the hxed point for the spin density wave critical 
found by Sur and Lee [36] using an expansion in e = 3 — d. They also also found the asymptotically 
nested Fermi surfaces in Eq. (1.3) under the 1-loop renormalization group flow of the e expansion. 
We will review their RG analysis in Section II. We then proceed to a computation of the optical 
conductivity in Section III, and find that the hot spot contribution obeys the hyperscaling of 
Eq. (1.1) (up to logarithmic corrections) in the e expansion, in agreement with Hartnoll et al. [31]. 
We turn to a computation of the non-zero temperature free energy density in the e expansion in 
Section IV. We hnd a result for the hot spot contribution to the specihc heat again in agreement 
with the hyperscaling of Eq. (1.1), and for reasons similar to those for the optical conductivity. 

Sections III and IV also examine the optical conductivity and the free energy in the limit 
of a vanishing bare 0 velocity: c —)■ 0. As the bare velocity is generically hnite, such a limit 
can only apply to observable properties over intermediate u or T: we hnd the allowed range is 
cA < u,T < vpA, where vp is a Fermi velocity (see Eq. (2.2)), and A is high momentum cutoff. 
Only in such a limit do we hnd hyperscaling violation as described by Eq. (1.2) with d* = 1. The 
quantum critical optical conductivity studied in Refs. 34 and 35 is analogous to this intermediate 
regime, and we maintain that their results do not apply when the the bare velocity c is not small. 

The more subtle question of the DC conductivity is examined in Section V; in discussions of 
the DC conductivity, we implicitly assume that uj T. Here, we have to consider the interplay 
between the hot spots on the Fermi surface with the remainder of the ‘cold’ Fermi surface more 
carefully [11, 31, 37, 38]. The cold fermions can short-circuit electronic transport, and so possibly 
dominate the DC conductivity. More generally, this belongs to a class of ehects associated with 
the conservation of total momentum, which can relax only via quenched disorder or umklapp 
scattering beyond that already continued in the continuum theory [11]. A general framework 
for describing such effects was presented in Refs. 5 and 14, using solvable holographic models, 
relativistic hydrodynamics, and memory functions. In the context of strange metals, it useful to 
begin with a microscopic model in which total momentum is exactly conserved [7, 11]. Then the 
conductivity can be written as [5] 




Q' 


A4 (—^cj)’ 


( 1 . 4 ) 
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where ctq is a finite and T-dependent ‘quantum critical’ conductivity, and the second term can 
be viewed as the contribution of the cold Fermi surface. The pole at cn = 0 has a co-efficient 
determined by static thermodynamic susceptibilities associated with the electric current J and the 
momentum density P, with Q = xjp and Ai = Xpp- These thermodynamic susceptibilities are 
usually non-critical, and so can be taken to be non-universal and T-independent constants, which 
depend on the full short-distance structure of the theory. Now we add perturbations associated 
with umklapp scattering or quenched disorder which can relax the total momentum [5-14, 39-43]: 
this leads to a momentum relaxation rate T which shifts the pole in Eq. (1.4) off the real axis to 
u = —iT, and so the conductivity takes the finite value at a; = 0 




Q2 1 

M (-ia; + r)' 


( 1 , 5 ) 


Note that T does have a singular T dependence associated with universal properties of the quantum- 
critical theory, and can be computed via memory functions [5-7, 11, 39, 44-47]. A notable feature 
[48] of Eq. (1.5) is that the quantum-critical (Tq and the momentum-mode conductivity are ad¬ 
ditive; this is in contrast to the Matthiessen’s Rule for quasiparticle theories, in which different 
quasiparticle scattering mechanisms are additive in the resistivity. The T dependence of the 
momentum-mode term in Eq. (1.5) was discussed in Ref. 11, using the assumption that the cold 
regions of the Fermi surface are ‘lukewarm’ i. e. the electron-electron scattering rate on the entire 
Fermi surface is faster than the impurity scattering rate; the results of Ref. 11 are not modified by 
the analysis of the present paper. 

Section V will present a computation of the quantum-critical conductivity ag for the case of a 
spin density wave quantum critical point in a metal in d = 2. The momentum mode contribution 
in Eq. (1.5) was computed in a previous work by two of us [11], and will not be addressed here. The 
computation of aq here is aided by the fact that the theory describing the hot spots is particle-hole 
symmetric. This implies that the scaling limit theory has Q = 0, and so we can cleanly separate 
away the momentum mode contribution; the full theory ultimately has Q ^ 0, but this arises from 
portions of the Fermi surface away from the hot spots [11]. Such a separation between aq and the 
momentum mode is more complicated in general [12]: in particular, for the Ising-nematic critical 
point there is no particle-hole symmetry to aid us, and we are not aware of any computation of 
aq for this case. For the spin density wave critical point, we compute aq in Section V using a 
Boltzmann equation method developed for conformal field theories [49-52]. We will carry out the 
Boltzmann analysis directly in d = 2, rather than the technically more cumbersome e expansion. 
Consequently, our results for aq will be qualitative, and not systematic. From the computations in 
Section V, we estimate that the leading T-dependence of aq has the same form as the cn-dependence 
of the optical conductivity: i.e. with bare velocities finite, hyperscaling is preserved with aq ~ 
constant; and with vanishing bare velocities, there is violation of hyperscaling with aq J'{dt-2)/z 
and dt = 1 and over intermediate T range cA < T < vpA. 
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II. MODEL 


In this section, we first recapitulate the low-energy continuum quantum field theory for fermions 
moving in a two-dimensional square lattice close to the transition to the antiferromagnetic phase 
with commensurate ordering wave vector Qaf = (tt, tt) [32, 33, 36]. We then explain the embedding 
by Sur and Lee [36] of the two-dimensional system into a higher-dimensional d = 3 — e space, and 
summarize the basic features of the e expansion. 

We begin by defining the action in frequency and momentum representation in two 

space dimensions x and y and one temporal (imaginary time) direction r: 




(m) , 


i=l m=± (T=t,i 

E 

1=1 



k J q 


(t>{q) ■ + h.c. 


(-)/ 

a' 


( 2 . 1 ) 


Here, the functional integral for the fermions goes over fermionic Grassmann fields ip, ip which 
carry additional labels according to their “home” hot spot (depicted in Fig. 1). Via a “Yukawa” 
coupling gf, the fermions are (strongly) coupled to a bosonic vector field with three components 
(j) whose fluctuations represent spin waves. At zero temperature kr is a continuous (imaginary) 
frequency variable with k = {kr, k) = {kr, k^, ky) and likewise for q. 

According to Fig. 1, the dispersions of the fermions e^(k) = ■ k in the hot regions are 


e]^(k) = -e^(k) = vp {vk^, ± ky) 

e^(k) = -ej(k) = vp (=Ffca: + vky) , (2.2) 

with V being the ratio of the velocities in x and ^/-direction; we will henceforth set vp = 1. In 
particular, the limit n —)■ 0 corresponds to locally nested pairs of hot spots, in which the Fermi 
line becomes orthogonal to the antiferromagnetic ordering vector Qaf and the fermion becomes 
one-dimensional and disperses parallel to Qaf- 

The physics of the action Eq. (2.1) in two space dimensions has been addressed with a vari¬ 
ety of techniques including resummation of subclasses of Feynman diagrams [32], field-theoretic 
renormalization group techniques [31, 33], and Polchinski-Wetterich flow equations for the effective 
action [53]. The bottom line is that the fermions and spin-waves are strongly coupled, one has to 
account for strong renormalization of the shape of the Fermi surface [33]. 

Here we embed the fermionic system in two space dimensions described by Eq. (2.1) into a 
higher-dimensional space; the “extra dimensions” are added perpendicular to the physical Fermi 
surface [36] that lies in the x-y plane and has co-dimension 1. Artificially introduced Fermi surfaces 
with co-dimension > 1 are gapped out by assuming a p-wave charge density wave order in directions 
perpendicular to the physical Fermi surface. This results in line nodes of the fermionic dispersion 
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with co-dimension 1 as needed. The main advantage of this embedding is that the density of states 
at the Fermi line is suppressed to p{E) ~ that is, it vanishes with energy for d > 2. This 

allows the powerful dimensional regularization techniques of relativistic systems to be adapted to 
the present problem. 

The d + 1-dimensional action 


4 Vc Nf 


j{k) [tr-K + + \ [ W + Tr [$(-g)<F(g)] 

n=l 17=1 j=l 

[3—d )/2 ^ r r 

+ ^^^7=—V y] y] / / [^n,<7,j{k + q)^„^„i{q)i-id-i^n,<7',j{k) + h.c] (2 

V^/ n=la,a' = li = l'^kJq 


is integrated over k = (K, k), which contains the physical momentum k = [kx, ky) and a d—1 = 2—e 
dimensional “generalized frequency” vector K = (kr,^) = {kr, ki,kd- 2 ), that includes the 
physical frequency kr in its first component and the d — 2 extra dimensions in the others and 
likewise for q. The bosons have been promoted to matrix helds (j){q) = '^^= 1 ^ (j)°‘{q)T^ with 
Tr[r“r^] = 25“^ conventions for the trace over S\J{Nc) generators r“. The fermions are collected 
in a SU(iVj-) flavor group and the physical limit of Eq. (2.3) is 


K —y kj- , 6 —y 1 , d —y 2 , — 2 , Nf — 1 . 


(2.4) 


Computations with Eq. (2.3) involve traces over products of d — 1 dimensional gamma matrices, 
collected in the vector (r, 7 rf_i) with F = ( 70 , f) = ( 70 , 71 ,..., 7 ^- 2 ), that satisfy {7;„7i.} = 2I6f,y 
and Tr J = 2. The book-keeping indices for the hot spots are: 1 = 3, 2 = 4, 3 = 1, 4 = 2; the 
two-component fermion spinors appearing in Eq. (2.3) disperse according to 

£i(k) = e+(k), 

£ 2 (k) = e+(k), 

£3(k) = e^(k), 

e 4 (k) = e^(k), (2.5) 


with the right-hand-sides defined in Eq. (2.2). The two-component spinors of Eq. (2.3) contain 
two of the original fermions from opposing sides of the Fermi surface [36]. 

Sur and Lee [36] performed a held-theoretic one-loop renormalization group analysis of Eq. (2.3) 
in d = 3 — e dimensions. They retained the simplest set of 5 independent running couplings. For 
the fermion propagator 3 wave-function renormalization factors are used, one in the direction of 
“time and extra dimensions” K and one each in the kx and ky directions. For the Bose propagator 
there are 2 wave-function renormalization factors, one in the Q direction and one for the qx,y 
directions (which have to be equivalent by point group symmetry). 
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The fixed point of the e expansion is defined in terms of the ratios A 


A ^ A* 
w ^ w* 


N^ + N^Nf-1 
N^ + N,Nf-3^ 

NcNf 


The fixed point determines a dynamic critical exponent 2 ; via 


g'^/v and w = v/c: 


( 2 . 6 ) 


2 = 1 + T + o(£2), (2.7) 

Note that and K scale as so the extra d — 3 spatial dimensions and the time dimension both 
scale as 2 ; with respect to the two physical dimensions, instead of just the time dimension as is 
usually the case with other models. 

While the scaling structure described so far is conventional, there are logarithmic corrections 
which arise from the flow of the velocities v and c flow to zero at long length scales. This flow is 
described by 


dc 

dln/i 

dv 


Az 


TT 


{z - l)c^ 
dc 


w 


d In /i d In /i ’ ^ ^ 

where /i is the renormalization group momentum scale. Such a dynamic nesting with n —)■ 0 
was found also in an earlier 1/Nf expansion [33]. At the fixed point with vanishing v and c, the 
antiferromagnetic ordering vectors intersect the Fermi surface at a right angle. This is illustrated 
in Fig. 1(b). Note that with n —)■ 0 at the fixed point, we must also have ^ 0 for the coupling 
A to remain finite; this is indeed found to be the case in the renormalization group flow. 


III. OPTICAL CONDUCTIVITY a{ui) 

In this section, we compute the optical conductivity (^{u) for fermions near the hot spots at the 
e expansion fixed point described in Section II. Our computation will be to order e, which requires 
evaluation of two-loop Feynman graphs. 

Before embarking on the description of the Feynman graphs, let us review the expectations of a 
general scaling analysis. The spatial directions, x, y, have scaling dimension 1, the time direction 
has scaling dimension 2 ;, and the 1 — e extra spatial directions with Dirac dispersion also have 
scaling dimension z. So the scaling dimension of the free energy density, F, is 

[F] = 2 + (2 - e)z. (3.1) 

The vector potential. A, has dimension 1, and so the electric current, J, being proportional to 
5F/5A has dimension 


[t/j — 1 -|- (2 — e)z. 


(3.2) 







Finally, the conductivity is given by the Kubo formula in terms of a current correlator, from which 
we deduce 

[a] = (1 - €)z. (3.3) 

These are the scaling expectations for a theory that obeys hyperscaling. If we have a violation 
of hyperscaling, we expect the spatial direction along the Fermi surface to not contribute in the 
counting of scaling dimension. So we should have 

[F] = 1 + (2 — e)z , [cr] = —1 + (1 — e)z, with hyperscaling violation. (3.4) 

We already know that the Fermi liquid contribution of the quasiparticles far from the hot spots 
violates hyperscaling as in Eq. (3.4) with z = 1. The question before us is whether the hot 
spot contribution preserves hyperscaling as in Eqs. (3.1) and (3.3), or violates hyperscaling as in 
Eq. (3.4). 

We hrst compute the one-loop (free fermion) contribution to the two-point correlator of the 
current density (JyJy). Then, we compute the two-loop (interaction) contributions to {JyJy), 
of which there are two: the “self-energy correction” (Section IIIB) and the “vertex correction” 
(Section IIIC). Finally, in Section HID we compile the results from the evaluation of the loop 
diagrams and, applying the Kubo formula, we derive the scaling form of the optical conductivity 
(t{u) for the fermions near the hot spots. 


A. One-loop contribution to {JyJy 


We have, for the current density in the y direction, 

j, = + jf + jf > + = 

Afc Nc 


cr=l j=l 


cr=l j = l 


(3.5) 


and likewise for J^ but with (1,3) -H- (4, 2). 

The one-loop contribution to this correlator is simply the non-interacting “bubble” containing 
a convolution of two fermion propagators as shown in Fig. 2, both from the same hot spot (we 
follow the index convention of Eq. (2.5) and absorb the identical contributions from the other hot 
spots into a prefactor) 


{JyJy)i-ioop{oo) = -2{l+v^)N,Nf I ^^^_^Tr[z7,_iGi(K,k)*7,_iG'i(K + W,k)], 


(27r)2 (27r)2- 


(3.6) 
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FIG. 2. Feynman graphs for the current-current correlator up to two loops. The black and grey boxes are 
current vertices for the n and n hot spot pairs respectively. The wiggly lines are the boson propagators 
and the solid lines stand for fermion propagators, (a): One-loop contribution for free fermions, (b): Two- 
loop self-energy correction computed in Section IIIB. There is also a partner diagram with the boson on 
the lower fermion line, (c): Two-loop vertex correction computed in Section IIIC. 


where W = (a;, 0) and the fermion propagator is given by 


G„(K,k) = 


r • K -|- 7d-i£n(k) 


(3.7) 


^ K2+£„(k)2 • 

We evalnate Eq. (3.6) using Feynman parameters in App. A 1 and obtain in Eq. (A2) to leading 
order in e: 


(4Jy)i_ioop(cu) = -Vl 4- j dk\\ NcNfU^ " (l^) ‘ 

where k\\ is the component of k along the Fermi surface of £i(k) (note fey = k^ for n = 0). For 
comparison with the subsequent two-loop contribution, it is useful to write this as 

{JyJyh-.o.M = -{l + V^)j If (yI) . (3.9) 

where the variable of integration £3 is a co-ordinate orthogonal to the equal energy lines of £ 3 (k). 

We can evaluate the integral over £3 to yield a factor of A, a large-momentum cutoff, and then 
we conclude that cri_ioop(i^) ~ We now observe that this result agrees with hyperscaling 

violating scaling dimension in Eq. (3.4) for z = 1. This is just the expected result, because we are 
dealing with the contribution of free fermions, and there is no distinction yet between the hot-spot 
contribution, and the Fermi liquid contribution of quasiparticles far from the hot spot. 


B. Two-loop self-energy correction {JyJy)sE 


To investigate the impact of interactions on {JyJy), we hrst compute the two-loop self-energy 
correction depicted in Fig. 2(b). There are two diagrams here with identical contributions, whose 
sum gives 


( Jy Jy)sE(^) 


-4:{l+V^)Nf 


rpi 

(27r)2 (27r)2-^ 


*7d-iGi(K, k)Si(K, k)Gi(K, k)*7,_iGi(K + W, k) 


(3.10) 


10 










FIG. 3. Key one-loop elements appearing in the two-loop self-energy correction (a) and two-loop (current) 
vertex correction (b). 


Here we note that this expression contains three fermion propagators from the same hot spot pair 
1 and one, inside the one-loop self-energy S„(K,k), depicted in Fig. 3 (a) from its “partner” hot 
spot pair 3. We now compute Si(K,k) separately. After that, we substitute the result back into 
Eq. (3.10) and perform the remaining integrations over k and K. The effect of large momentum- 
transfer scattering of fermions from one hot spot pair (1) to its partner (3) via exchange of bosonic 
spin fluctuations is captured in the self-energy 

Si(K, k) = ^ ^ (rV^) J + K, q + k)i'yd-iD{Q, q) , (3.11) 

where the spin fluctuation propagator 

involves the spin-wave velocity c which vanishes at the Sur-Lee hxed point near the hot spots, as 
does the Yukawa coupling the ratio jc however attains a hnite value (see Sec. II). 

We evaluate the expressions hrst in Appendix A 2 using a simplifying approximation valid only 
for small bare velocities c and v. In the limit v,c ^ 0, the integrand in Eq. (3.11) then depends 
on Qx only via the spin fluctuation propagator. Thus, we can hrst perform the integration and 
then set c = 0, which is equivalent to replacing 

f V [ 

J ( 27 r )2 J ( 27 r )2 

^(Q,q)^-T^- (3.13) 

c IQI 

in Eq. (3.11). This way, Eq. (3.11) picks up the hnite prefactor g'^/c and the integrand becomes 
independent of both velocities v and c. The resulting integrals are performed in Appendix A 2 
using Feynman parameters to obtain 


Si(K,k) = Si(K) = 


TT^ ^/2r(e/2) g- 


-i- 


(2vr: 


4-e 


N. 


N^-l 

Y 




f C ^' 

J .7 = 1 


r K 

[ K 2]^/2 


Iq 


1_€ 
2 2 


(3.14) 
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We observe that the above fermion self-energy depends only on frequency and not an spatial 
momenta k anymore; the fermions near the hot spots see an essentially “local” boson with the 
“1 over frequency” propagator of Eq. (3.13). Eq. (3.14) induces an anomalous scaling for the K 
dependent part of the fermion propagator, but the absence of anomalous dimensions for the spatial 
k components renormalizes the dynamical exponent z to values larger than one at the Sur-Lee hxed 
point. 

However, for our purposes here, the approximation associated with Eq. (3.13) turns out not to 
be sufficient, since v and c vanish only logarithmically near the hot spot. It is thus crucial to obtain 
the full V and c dependence of the pole term in Eq. (3.14), and of that in Eq. (3.10). The needed 
integrals are computed in Appendix C, and the hnal result for the two-loop self-energy correction 
is the rather complicated expression in Eqs. (C2) and (C3). Computing its singular pole in e and 
dropping power divergent terms, we obtain 


{jyJy)sE{,^) ~ 


des (A^c-l)^Va^^ ' , {I - + v'^) / 2 \ 

2v Gdvr^ce Jq ^ (c^ -|- x(l -f- — c^))^/^ ^2 — c^) / 

(3.15) 


in terms of the same £3 variable of integration used in Eq. (3.8). 


C. Two-loop vertex correction to {JyJy)vevt 


The vertex correction graph. Fig. 2(c), is considerably more involved than the self-energy cor¬ 
rection of the preceding section; although for c —)■ 0 it is free of 1/e poles of the type Eq. (3.15). 
Using the abbreviation Si(K, k, W) for the one-loop current vertex correction in Fig. 3(b), we can 
write the entire graph including contributions from all hot spot pairs as 


(JyJ„)vert(w) ^ -2(1 - V^)iNf 


r 

f d^-^K 


/ (27r)2j 

' Tr 

(27r)2-^ 

7,_iGi(K, k)S3(K, k, W)G'i(K + W, k) 


(3.16) 


We observe that Eq. (3.16) contains two fermion propagators from one hot spot pair and two 
fermion propagators (inside the current vertex Si) from the partner hot spot pair, unlike the self 
energy correction Eq. (3.10). The one-loop correction to the Jy vertex. 


(K,k,W) = 


2 „. 




Nf — 

} j=i 


^(rV^) J 


(27r)2 (27r)2-^ 


7 <i-i^ 3 (K -|- Q, k -|- q) 7 rf_iG 3 (K -|- Q -|- W, k -|- q)7d-i pj 2 1 2 2 

“T C 

(3.17) 


.72 
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does not contain a 1/e pole in the limit of c —)■ 0. This subsequently leads to the lack of a pole 
in the two-loop vertex correction (the details of the computation are presented in Appendix A3). 
For n, c 7 ^ 0 however, as described in Appendix C, the vertex correction picks up a small pole with 
a coefficient of 0{c): 


(A/A/)vert (^) ~ 


desiN^^ (1 - 

2v 327lh Jo "" (c2 -h x{l +v‘^- c2))3/2 + x{l + - C^) 

(3.18) 


However, this is subdominant to the self energy correction in Section IIIB, as the latter is hnite 
in the limit c —)■ 0; this is a consequence of a Ward identity discussed in Appendix C. 


D. Renormalized conductivity cr(uj) 


We can now add the leading free contribution in Eq. (3.9), the singular self-energy correction 
in Eq. (3.15) and the appropriate counter term to obtain the renormalization of the conductivity 
(t{u) near the hxed point 


ayy{uj) ^ (1 + n^) 


dso Nf-Nf _ 

- -UJ 

2v Ibvr 



1 + 


{z-l) {l-xf/'^ (u? 

TT (c^-f x(l-I-— c^))3/2 ^ -I-x(l-f 



(3.19) 


The interpretation of this central result requires some care in the limit of small v and c, and we 
consider various cases separately below. The important point here is that the argument of the 
logarithm is of order (o;^ -|- (c^s)^)///^, and so the renormalization-group-improved perturbation 
expansion will lead to powers of (ca^ - 1 - (c^a)^)///^, in contrast to the power of oj alone outside the 
square bracket. This difference arises because, at leading order, the singular contribution of the 
hot spot is the same as the rest of the Fermi surface, while at higher orders there is quasiparticle 
breakdown only close to the hot spot. Consequently, there is a modihcation in the nature of the 
£3 integral, where we recall that £3 measures distance away from the hot spot along the Fermi 
surface. 

First, let us assume the bare value of c is so small that the £3 dependence of the argument of 
logarithm can be ignored; this was, effectively, the limit that was implicitly taken in by Abanov 
et al. [34]. This requires that cA < ca, where A is the momentum space cutoff. Then, we can 
easily perform the integral over x, and after re-exponentiating the logarithm in the e expansion, 
we conclude that 


de 




3 -e 
CJ 


A/i 


2v 


l-Z ^-€+(z-l)^ 


1 + {z -1) ln(a;//i) 


(3.20) 
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This is the answer expected from the hyperscaling violation case in Eq. (3.4) at this order in e; 
Abanov et al. [34] fonnd (y{u:) ~ 1/\/uj, which is consistent with Eq. (3.4) for their dynamic critical 
exponent z = 2. Thns the uj dependence of a violates hyperscaling as in Eq. (1.2) with dt = 1 for 
cA < u < A, which can be within a nniversal regime only if the bare valne of c is small enongh, as 
claimed in Section I. 

Next, we consider the more generic case where the bare valne of c is of order nnity. Then, we 
can divide the integration over £3 in Eq. (3.19) into two regimes. There is the far from the hot-spot 
regime where ce^ ^ ca, and the close to the hot-spot regime of ce^ ca. The contribntion of the 
close to the hot-spot regime is similar to that in Eq. (3.20), except that the npper bonnd on the 
integration over £3 involves u 


a. 


yy 




de. 


2v 




(3.21) 


Actnally, by scaling, we expect the npper bonnd on the momentnm integral over e^/v to scale as 
at higher order in e; nsing snch an npper bonnd in Eq. (3.21) we obtain the generic hot-spot 
contribntion 


^yy\^) 


UJ 


-e+(z-l) + l/z 


(3.22) 


To this order in e, this is the scaling expected by the hyperscaling preserving scaling dimension in 
Eq. (3.3). In comparison to the hyperscaling violating answer obtained in the direct n, c —)■ 0 limit 
nsed for Eq. (3.20), the condnctivity has acqnired an extra factor of So we reach one of onr 
main conclnsions, that the hot-spot contribntion to the condnctivity generically obeys hyperscaling 
as in Eq. (1.1). We have not written ont explicit factors of v and c in the final scaling forms, bnt 
these are nltimately only expected to yield powers of (ln(l/a;))“\ and so hyperscaling is only 
obeyed np to powers of ln(l/ci;). 

Finally, we also have to consider the contribntion of the far from the hot spot regime ce^ ca. 
In this regime, the term inside the sqnare brackets in Eq. (3.19) is ^-independent, and so we 
obtain an additional contribntion a ~ This is jnst the additive Fermi liqnid contribntion of 
long-lived qnasiparticles far from the hot spot. 


IV. T > 0 FREE ENERGY 

In order to stndy the finite temperatnre dynamics of this model, we need to compnte the free 
energy density at T > 0. The free energy density has contribntions from the free fermions, the free 
bosons, and a “self energy” correction dne to their interactions. Following the lessons learned in 
the analysis of the optical condnctivity in Section III, we will perform the compntation here only in 
the simpler limit of vanishing velocites n, c —?■ 0 , where we can replace the boson propagator by the 
momentnm-independent form in Eq. (3.13). However, as described in Section HID, we will assnme 
that the low T hot spot contribntion for the case of finite velocities can be estimated by limiting 
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the range of the fermionic kx integral (along the Fermi surface) to an upper limit we 

have assumed the the upper cutoff is determined by T rather than u for the optical conductivity 
in Section HID. 

The free fermion, Fj?, and free boson, contributions to the free energy density, F, are obtained 
straightforwardly to leading order in e (the prefactor of 4 in the fermion contribution comes from 
having 4 pairs of hot spots) 


F = AN^NjT ^ 


dkx f dkyd} , 

^ In 


2tt 


2-t 


= / dkx N,NfT 


3-e 


{2n) 

(Km 

V 27r2 ’ 


(1 + g(fe 2 +K 2 )l/ 2 /T^^^ ^ ^_(fc 2 +K 2 )l/ 2 /Ty 


(4.1) 


where the inhnite temperature-independent constant part was dropped. For the bosons 


F“ = (1 - ]V2)T 



[ In 

'l _ g-(c2q2 + Q=)l/Vr 

/ mJ 

(27r)i-^ 

- 


TT 


90c2 


(N‘‘ - 1)T 


4-t 


(4.2) 



FIG. 4. The simplest interaction contribution to the free energy at 0{g^ 


The interaction contribution to the free energy at two-loop order is given by Fig. 4. It may be 
expressed as 


Ffb = T^Tr 


V2-l 

i=i 






3-e 




n(«.r) 


Q2 -H c2|q|2 a;2 ’ 


(4.3) 


where n(g, T) is the RPA polarization bubble given by 

d^kd^-^K 


n{q,T)^g^J 2 

n ' 

We separate out n(g, T) as 




3-e 


T ^ Tr ['yd-iGn{k)'yd-iGn{k g)]. 


(4.4) 




n(g, T) = (n(g, T) - n(g, 0 )) + n(g, 0 ), (4.5) 

and evaluate the finite temperature part setting n = 0 at the outset, taking g,c ^ 0 with g^/c 
hnite and equal to its fixed point value X*w*. The zero temperature part is evaluated with v 0 
at the outset, and g, n —)■ 0 with g‘^/v finite and equal to its fixed point value A* . As described in 
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Appendix B, this separates out the contributions that renormalize the free fermionic and bosonic 
contributions, with the hnite temperature part of 11 renormalizing the free fermionic contribu¬ 
tion and the zero temperature part renormalizing the bosonic contribution. We then obtain (see 
Appendix B for details), for the singular parts. 


Ffh — -I- 

~ ^ fb + ^ fb 5 

Fjl'> = j dK {Nl - 



{Nl - 


360nc^e 


3C(3) \ 
Ibvr^e ) ’ 


(4.6) 


We have set the momentum renormalization scale /r = 1 in the present section. We thus get, after 
plugging in the hxed point values. 


Ff = F° + F^i^ = / dF 


fb 


3C(3) 

27r2 


1 + (^ - 1 ) 




2F 


(4.7) 


where the pure 1/e pole is cancelled by the usual addition of a counter term. Similarly, 


n = F° + f)? = - 1)7“- 


i + 2i^r- 


-)■ 


TT 


90c2 




We now observe that the bosonic term F^ is compatible with the behavior rsj T2 ^+2/^ expected 
from the hyperscaling preserving scaling dimension in Eq. (3.1); the agreement holds to hrst 
order in e after recalling that 2; — 1 is 0{e) from Eq. (2.7). For the fermionic contribution, as in 
Section HID the behavior depends upon the fate of the integral. As noted at the beginning 
of the present section, for the low T behavior we should impose an upper cutoff on the integral 
of order then Ff ~ 2"3-e-(2-i)+i/2 agrees with ~ ^ hrst order in e. 

Thus both the bosonic and fermionic contributions to the free energy obey hyperscaling, and the 
behavior in Eq. (1.1), up to logarithms. 

As was the case in Section III, for very small bare velocity c, and for cA < T < A, there is 
a regime of hyperscaling violation when the kx integral is replaced by A, and behavior is as in 
Eq. (1.2). Note that we are using units in which the velocity vp in Eq. (2.2) has been set equal to 
unity; so the full condition for this intermediate regime is cA < T < vpK. 


V. QUANTUM BOLTZMANN EQUATION 

We now compute the hot spot conductivity aq appearing in Eq. (1.4) in d = 2 using a quantum 
Boltzmann equation approach [49-52]. We use the Keldysh formalism at one-loop order to derive 
quantum kinetic equations for the fermions and bosons in the presence of an applied electric held, 
and then solve these equations in linear response to obtain the contribution of the fermions near the 
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hot spots to the DC conductivity. Note that, unlike the previous sections, we are not performing 
a systematic e expansion here, but working directly in d = 2 to minimize technical complexity. 

As in Section IV, we will restrict our analysis to the case of vanishing v and c, when the Fermi 
surfaces are nested, and manipulations similar to Eq. (3.13) can be applied. With hnite v and c, 
as argued in Sections HID and IV, we can estimate the low T hot spot conductivity by limiting 
the kx integral along the Fermi surface by an upper bound of order 


A. Keldysh framework 

We begin by expressing the action in Eq. (2.1) on the closed time Keldysh contour [51, 52]. 
Denoting with subscripts + the forward part of the contour and with subscripts — the backward 
part of the contour, we obtain for the free part of the action 


A, 


'Ipip 




(2ir) 


2/^*7 (2^)2 


£=1 m=± (T=tq 


p) - er(p)) p) - p) - er(p)) p) 

-q) • {-dt - Wq) q) - -q) • {-dt - ^q) q) 


(5.1) 


with ojc^ = c|q|. The interacting part is given by 




/ OO ^ 4 

dt / dhj2 

■°° i=l t7,cr'=t, 


r) • r)- 


(-) 


£=1 o-,o-'=t,4. L 

d+) j+ 


- r) ■ i/'v:-(i,x)0(i,r)i/iy __(i,r) + h.c. 


(5.2) 


We now perform the standard bosonic and fermionic Keldysh rotations: For the real bosons we 
use 


4>+ — 4>c + <Pq, 


and for the Grassmannian fermions we have 


~ ^We,cr,l ^ V£,o-, 2 h 


, (m) _ I (m) _ , M X 

k’e,a- ~ ^We,cr,i Fq(T,2h 


(11 

c ~ 




~ ^W£,cr,l W iP£^cr,2)l 
- 'I 


(5.3) 


(5.4) 
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Hence, we get for the free fermion part of the lagrangian 


I -1 


i=\ m=± (T=t,4, 


{Gf^\ 

0 




V'Sl(^)P) 


-1 


(5.5) 


where the inhnitesimal 5j ensures convergence. Inverting this matrix, we obtain the bare fermion 
Green’s function matrix 


G £m _ 

0 “ 


G BIm ^Klm 
0 ^0 

0 


(5.6) 


The bare retarded (R) and advanced (A) fermion Green’s functions thus are 

1 


Gf™(a;,p) = 


Gf”^(a;,p) = 


cj + f0+ - e7*(p)’ 
1 


cj - f0+ - e^(p) ‘ 
For the free boson part of the lagrangian we have 

0 

[D§r '5f 


(5.7) 


-1 


= 2 i “h) -q) 


0c(w,q) 

q) 


(5.8) 


where the inhnitesimal 5^ again ensures convergence. After performing the matrix inverse. 


Dn = 


Di 0 


(5.9) 


and the retarded and advanced boson Greens’ functions hence are 


q) = ^ 


D^{u, q) = 


2 (cj + iQ+f — 
1 

{uj - i^+f - 


(5.10) 


The interaction between fermions at the (£, +) and (£, —) hot spots and the boson takes the 
following form; 




"'(p'lp'lp 




(+) 


0c(i,r)-W 0g(t,r)-w \ / V’y,i(f,r) 


(j)q{t, r) • w r) • I \ r) 


i=l (t,(t'=|,4. 


+ h.c. 
(5.11) 


This gives rise to the Feynman rules summarized graphically in Fig. 5. 
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FIG. 5. Feynman rules and one-loop graphs for the self energies of the fermions and bosons in the Keldysh 
formalism: (a) Fermion propagators, (b) Boson propagators, (c) Yukawa vertices, (d) Self energies. Here 
X = (t, r), hot spot indices (^, ±) are suppressed and the external legs on the self energy diagrams are 
amputated. The external momentum and frequency on the self energy diagrams is on shell. The 2-1 and 
q-q propagators are zero and hence are omitted. 

We adopt the shorthand convention oix = (t, r) and q = (m, q) to combine spatial and temporal 
components. We have the relations 

QKim ^ Qmra ^ ^ QAlm ^ 

= D^oFb-FbO ( 5 . 12 ) 

where C = Ao B implies C{x,x') = f dxiA(x, Xi)B(xi, x') and Ff^b are respectively the fermionic 
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and bosonic distribution functions. The Dyson equations for the matrix fermion and boson Green’s 
functions are 

o = 1, 

-Eb)oD = l. (5.13) 

The self energy matrices S have the same form as the inverse Green’s function matrices in Eq. (5.6) 
and Eq. (5.8), and hence the different components of the self energies are given by the graphs 
in Fig. 5 at the one loop level. Dehning central and relative coordinates Xc = (x + x')/2 and 
Xr = {x — x')/2, we can convert the two point functions G, D, F and S which are of the form 
A{x^x') = A{xc + Xr/2^Xc — Xr/2) iuto a momentum representation via the Wigner transform 

A{xc,p) = j dxre~"^''^A (xc + y,a:c - yj . (5.14) 

Since we have spatial translational invariance in the linear response limit a of weak applied electric 
held E, we can further simplify A{xc,p) A{t,p). We will also always consider external particles 
to be on shell in the subsequent computations of the collision integrals. We dehne an alternate 
parameterization ff^j, of the distribution functions Ff^b 

Ff{t,p,uj) = 1 - 2ff{t,p,u), 

Fb(f,q,n;) = 1 + 2/b(t,q,a;). (5.15) 

In thermal equlibrium in the absence of any applied electric helds, we have //_6(f,k, w) = 
where = 1/(1 ± are the thermal Fermi and Bose functions respectively [52]. 

B. Kinetic equations for fermions and bosons 

There are two coupled quantum kinetic equations [52, 54], one for the electrically charged 
fermions, 

(^ + E. T) Fp{t,p) = (6.16) 

with the on shell fermion distribution function Fj^(t, p) = Ff(t, p, e/"(p)); and one for the neutral 
bosons, 

yF,(f,q,a;q) = /.“"[Ey, E,](f, q), (5.17) 

The fermion electric charge is set to 1 for simplicity. The two collision integrals have the general 
form [52] 

m, P) = Pi e'^(p)) + 2Ef (f, p)Im [Sf ±(t, p, e'±(p))] , 

=fSf(f,q,a;q) + 2F6(f,q,a;q)Im[Sf(f,q,a;q)] . (5.18) 
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At the one-loop level, are given by the graphs in Fig. 5. The self energies and collision 

integrals are computed in Appendix D. We obtain 


iffl[ff,h]{t,p) = 


= 3g^ 


d\ 1 
271 Un 


<5(ef(p) -ej{p-q) - u^) p)(l -p - q)) - P - q)/b(t, q, a;q) + 






+ ff'^{t,p)fb{t,q,uj^) -6{ef{p)-eJ{p-q)+uj^) p)(l -/^ (t,p-q))- 




e±i 


p - q)/fe(t, q, -Wq) + p)fb{t, q, -Wq) 


'e±/ 


(5.19) 


and 


4"‘i//.A]((.q) = 49^E / 

i 

+ k + q)fb{t, q, Wq) - 


d^kr 


5 {cf, (k) 


-e/(k + q)) (^/^+(t,k + q)(l -/^ (^,k)) + 


// k)/b(t, q, Wq) j-I-(-1--H-—) , 


(5.20) 


where we have expressed the collision integrals in the alternate parameterization (5.15) of the 
distribution functions Ff^b and fj^{t,p) = //(t, p, e^(p)). 


C. Ansatz and solution for conductivity 


If we set the collision integrals to zero, the distribution function for the neutral bosons unaffected 
by the applied electric held remains hxed at its equilibrium value. For the Fermions, we have 

+ (5.21) 

To solve this in linear response, we switch from the time to the frequency domain and parameterize 
the deviation of /j^(a;,k) from its equilibrium value by [54, 55] 

p) = 27rh(a;)ny(v± ■ p) + ■ E(a;)(/?(v± • p,a;)n/(v^^ ■ p)(l - n/(v± • p)). (5.22) 


Inserting this parameterization into Eq. (5.21), we obtain the collisionless Lf function in linear 
response 


^nc 



■p,a;) 


1/T 

—ioo 0 +' 


(5.23) 


We have the electrical current density: 


JM = 2$^ [ -|j^vr(V’rVr)(^iP) = 2$^ 5^ / 

l=\ m=± ^ ' e=l m=± 


d^p 


vr/^(a;,p). 


(5.24) 
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and hence obtain the linear response conductivity 



2EE / 

n —1 ™—u ^ 


4 



£=1 m =± 


It is easily seen that if the 9?nc is used in the above expression, we obtain a temperature 
independent collisionless conductivity that has a delta function in u in its real part. In fact, for 
the Sur-Lee embedding in higher dimensions, we have for the collisionless conductivity 



(5.26) 


which is derived in Appendix D; assuming the kx integral yields a factor of the cutoff A, this yields 
a conductivity with the hyperscaling violating scaling dimension in Eq. (3.4), as expected for free 
fermions. Once collisions of the fermions with the bosons are included, these delta functions are 
broadened, and the k^ integral has to be performed with more care, as in Section III D. 

Returning to d = 2, we hnd that, to linear order in E, the bosons still remain in equilibrium 
and their distribution function is hence given by the thermal Bose function /^(t, k, Wk) = '^■fe(to'k) if 
the parameterization Eq. (5.22) is used (See Appendix D for a derivation of this fact). Intuitively, 
this is because the linearly dispersing hot spot model exhibits particle-hole symmetry, making the 
charge carrying modes excited by the applied electric held particle-hole pairs with the particle 
and hole moving in opposite directions. The bosons then do not absorb any momentum that they 
have to dissipate when the particle-hole pairs recombine and hence remain in equilibrium. This 
behavior is also present for quantum critical transport in graphene [54]. 

We now consider the system to be at the hxed point discussed previously, in the spirit of [49]. 
We take the applied electric held to be in the y direction (E = Eyey); since n —)■ 0 at the hxed 
point, only the i = 1 and £ = 3 pairs of hot spots contribute signihcant response in this case. (For 
the electric held in the x direction we obtain the same response with the £ = 4 and i = 2 hot 
spot pairs respectively). We insert the functions parameterized by Eq. (5.22) and the thermal 
Bose function for fh into the frequency domain version of the fermion kinetic equation Eq. (5.16), 
and linearize in Ey to get the following integral equation for ip in the n —)■ 0 limit (considering the 
£ = 1 pair of hot spots) 
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13 c /2 r (Pa 


5{2py -Qy- c|q|) p{py,uj)nf{py) (l - (l - Ti f {Qy - Py) ^ nfe(c|q|)) + 


2 271 J c|q| 

+ ^i%-Py,^)nf{qy-Py) {1 - Ufiqy - Py)) {nf{py) + nb{c\q\)) 

- <5(2py -qy + c|q|) p{py,uj)nf{py) (1 - nf{py)) (1 - - py) + nb(-c|q|)) + 

+ Piqy-Py,^)nf{qy-Py) {I - nf{qy - Py)) {nf{py)+nb{-c\q\)) 


= {-iU + 0+)p{Py,Uj)nf{Py) (1 - UfiPy)) - -HfiPy) (l - H f (Py)) 


(5.27) 


In the collision term, the boson momentum parallel to the Fermi surface, qx, is limited by the Bose 
function to a value of order T/c. Integrating out qx, we obtain 

^9^ 


2770 


CoiPy)(p{Py,Uj) + 


Cl[^,Py] 

^fiPy) i^-nf{Py)) 


= {-iuj + OpcfiPy, Uj) - (5.28) 


where. 


1 

CoiPy) = 2 / dqy 


sgn(g^ - 2py) — (1 - npqy - Py) + nyippy - qy)) 


^{qy - 2py)^ - c^q, 


Ci[p,Py\ = 


dqy 


sgn(gy - 2 py )— ip{qy-py,uj)nf[qy-py)x 


\Ji.qy - ‘^Pyf - c\ 


X (1 - nf{qy - Py)) {nf{py) + UbiPpy - qy)) 


(5.29) 


This equation may be solved iteratively by choosing the trial solution 


ipl{Py,Uj) 


-1/T 

^^CoiPy) + iu - 0+ 


(5.30) 


and iterating 


2710 


CoiPy) +iuj - 0+^ ipj+i{py,u) 


1 _ 3p^ Ci[(pj,Py] 

T 27ionf{py){l-nf{py))' 


(5.31) 


for j > 1 (Note that pi in Eq. (5.30) may be derived from inserting (po = 0 into Eq. (5.31)). The 
integral for Ci may be evaluated numerically by sampling (p at a discrete set of points and then 
constructing an interpolating function through these points. The exponential decay of nj^b at large 
values of their arguments ensures convergence of the integrals and supresses errors arising from the 
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extrapolation of the interpolating fnnction to large arguments. The trial solution (/?i is actually 
fairly accurate, and this algorithm converges in a small number of iterations. 

In the limit of c —?■ 0, which also occurs at the hxed point, we see that Cq ~ 1/c because the 
singularity in nb{2py — qy) as 2py — —)■ 0 is cut off by c in the 0 function in Eq. (5.28); this also 
occurs in the integral for Ci in (5.29). Hence in this limit, we have for a; ^ 0, 

X'(y) (5.32) 

for some function y. It can be seen from Eq. (5.31), and also established numerically, that given 
the above form for (pi, all <^j>i will also be of the same form. Numerically, we hnd that y is an 
even function and y(0) = 0. Since y is an even function, it is easy to show that the (1, —) hot 
spot contributes the same value to the conductivity as the (1, +) hot spot as it is related by the 
transformation of Py —?■ —py in the above computation. Similarly, the (3, ±) hot spots also produce 
identical contributions equal to those from the (1, ±) hot spots. Hence, using Eq. (5.25), 

y (^)^/(Pj/)(l-’^/(Py))) ~ (5.33) 

In the last step we have changed the fermion momentum notation from px to kx for compatibility 
with earlier discussions. It is also easily seen that (Jxx = (^yy if we repeat the above analysis for the 
£ = 2 and £ = 4 hot spots instead, and that axy = CFyx = 0. So Eq. (5.33) is the estimate by the 
Boltzmann equation of the value of the conductivity (Tq in Eq. (1.4). 

We now need to determine the T dependence implied by Eq. (5.33) using scaling ideas. Under 
the renormalization group flow, we expect that the coupling A = /v flows to a hxed point value. 
While this hxed point can be determined precisely under an e expansion, we are only able to make 
an estimate in the present computation carried out directly in d = 2, where g^ jv is a dimensionful 
quantity of order . The natural scale for the momentum p is set by the temperature, and so 
p ~ So in d = 2, we can expect that g'^/v ~ Ignoring the logarithmic factors, we 

therefore have the estimate 

ag ~ y dkx T-Vt (5.34) 

Finally, as in Sections HI and IV, we bound the kx integral by to conclude that aq ~ constant, 
as claimed in Section I. And also as in previous sections, for a small bare c, we will have aq ~ 
in the intermediate T regime cA < T < A (and as noted at the end of Section IV, after 
restoring units, this condition is cA < T < vpA). 


VI. CONCLUSIONS 

We have computed the critical conductivity and free energy at the onset of spin density wave 
order in metals in d = 2 using the e expansion introduced by Sur and Lee [36]. The advantage 
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of this method is that the e expansion appears to be valid systematically order-by-order in e, and 
there is no breakdown in the renormalization group flows. The e expansion exhibits a logarithmic 
flow of the velocity v to zero at large length scales, and a dynamic nesting of the Fermi surfaces. 
We found that hyperscaling was obeyed, with the hot spot contributions scaling as in Eq. (1.1). 

It is interesting to compare these results with a previous two-loop, large N renormalization 
group analysis of the spin-density wave critical point in Ref. 33, which also found a logarithmic 
flow of V to zero at low energies. However, it was also found that the large N expansion broke 
down at sufficiently large scales. The same large N framework was used to compute the optical 
conductivity by Hartnoll et al. [31], and it was found that hot-spot contribution was (t{u) ~ 
constant in the limit of vanishing n, as expected under hyperscaling in d = 2. 

This similarity between the large N and e expansion indicates that the terminology ‘quasi-local’ 
for the latter expansion [36] should be used with some care, and we have avoided it here. The 
basic scaling properties are similar to those of a standard, spatially-isotropic, critical point obeying 
hyperscaling with a hnite dynamic critical exponent 2; given by Eq. (2.7). The deviations from 
strong scaling arise only in logarithmic corrections, which are linked, ultimately, to the asymptotic 
nesting of the Fermi surfaces [32, 33] in Fig. lb. 

We also carried out computations for the free energy density at non-zero T using the e expansion. 
Again our results were in excellent accord with hyperscaling expectations. Both the fermionic 
excitations at the hot spot and the collective bosonic fluctuations scaled with the same power of 
T, as shown in Section IV. 

There was, however, for the somewhat unnatural case of a sufhciently small bare boson velocity, 
an intermediate energy regime where hyperscaling was violated. This was discussed in Section III 
for the optical conductivity, and in Section IV for the free energy. The optical conductivity results 
of Refs. [34, 35] are similar to this hyperscaling violating regime, and our analysis indicates that 
their results do not apply when the bare boson velocity is not small. 

Finally, in Section V, we addressed the question of the DC conductivity. Because of the con¬ 
straints of total momentum conservation, such a computation must be carried out [5, 7] in the 
context of the expression in Eq. (1.4), which separates a quantum critical conductivity ctq from 
that associated with ‘drag’ from the conserved momentum. We estimated (Tq in Section V and 
found a result that scaled as up to logarithms. Thus, the aq contribution to Eq. (1.5), in the 
theory of the spin density wave critical point, is likely not the mechanism of the strange metal 
linear resistivity. 

The momentum-drag term in Eq. (1.5) was considered in a previous work by two of us [11] for 
the spin density wave critical point: there we found that quenched disorder which changes the 
local critical coupling did lead to a linear-in-T resistivity. This conclusion is not modihed by the 
considerations of the present paper. 
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Appendix A: Computation of (JyJy) 


All computations in this appendix are symbolically described for the n = 1 fermions, with the 
identical contributions for n = 3 accounted for by doubling the overall prefactors. Momentum 
integrals in dimensional regularization are performed using the standard identity 


f M,. ^ r((a + d)l2)T{b - (g + d)/2) 

J (P + A)* V{dl2) r(6) 


1. Free fermion contribution to {JyJy) 


The free fermion contribution to {JyJy) is given by Fig. 2(a) and is straightforwardly computed 
in dimensional regularization: 


( Jy Jy ) free (^ ) 

- 2(1 + v‘^)NcNf 


at I 


(27r)2 (27r) 


2-e 


■. , ..r-K + 7,_i£i(k). , .,r-(K + w) + 7rf_i£i(k) 

+ (K + W)^+ei(k)^ 

-K- (K + W) + £i(k)2 


= -4iViv,vTT^ fdk, Cdx _ 

c /V +v J a ax J (27r)2-[(K + xW)2 + £i(k)2 + a:(l-x)W2]2 


d^-^K f -K‘^ + W^x{l - x) 


^ NcNf , f T - I ~~ ' •• 

27r y "70 7 (27r)2-^ V[K2 + a:(l-a:)W2]3/2 ^ [K2 + x(1-x)W2]V2 


NcNfir 


l-e/2 


7rF(l — e/2)(27r) 


2-e 


[ dk\\ [ dx(^ -- ^{x{l-x)Y/^-^/^+ 


'o Y 


2F(l/2) 


F(1 6/2)F(1/2 + 6/2) _ ^^^l/2_e/2 _ F(2 e/2)F( 1/2 + 6/2) ^^^^ _ ^^^l/2-e/2 

2r^3/2^ 2r^3/2^ 


= -x/TT^ j dfc|| + 0{e)^ , 


(A2) 


where k\\ is the component of k along the Fermi surface. 
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2. Fermion self-energy correction to (JyJy) 


In this subsection, we will freely take the limit of vanishing velocities associated with Eq. (3.13). 
The extension to the case of hnite velocities will be presented in Appendix C. 

This self-energy correction is given by Fig. 2(b) and a partner diagram with the boson on the 
lower fermion line. The sum of the two gives 


(JyJy) SEi^) 

= 4iNf 


d^k d^-^K 


(27r)2 (27r)2- 


Tr 


7d-i 


r ■ K + ^d_iky r ■ K + 'jd-iky r ■ (k + w) + ^d-iky 


K2 + fc2 


-Si(K,k): 


K2 + A;2 


-Id-r 


(K + W)2 + kl 

(A3) 


We hrst compute the fermion self energy, given by Fig. 3(a): 


Si(K,k) 

/u2_i 

= A: w/dATrA( r.(K + Q)-T.-.fa + fc,) _ 1 

N, J ' (K + Q)2 + + 

^2_1 

= _ -r-K(l -x) - 7 rf_i(g^ + fc^) _ 

Nf ^ Jq J (27r)2 (27r)2-^ [Q2 + x(l - x)K 2c2q2(l - x) + a:(gj^ + fcy)2]2 

pj2 _-I 

_ 7r^-^/^r(l + e/2) f j j. f\ f _ -T ■K{l-x)--fd-i{qy + ky) 

(27r)2-'= Nf Jo J (27r)2 [a:(l - a;)K 2 4- c 2 q 2 (l - x) -1- a;(gj^ +/Cy) 2 ]i+A 2 

j\^2 _1 

„3/2-,/2r(i/2 + £/2)jV ^7, , „ fy f dq, -T-K{l-x)-'u.dq, + K) 

(A4) 


where in the last step we integrated out and then sent c —)■ 0. After shifting Qy by ky we integrate 
it out to get 


Si(K,k) = -* 


„2-,/2p(j/2) g2^, 

124 


Af?-1 


4-6 


cNf — 

J j=l 


(rV-^) / dx{x{l — x)) 


-72 


1-x r-K 


X 


[K2]72' 


(A5) 
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Inserting this into the expression for {JyJy)sE-i 


{'4'4)se(w) 

= 16(1 - N'i) 


^^^ 5/2-./2r(,/2)r(l/2-e/2)^V 

2i-^(27r)8-2T(l - e/2) c 


= 32(1 - Nl) 


7r5/2-^/2r(e/2)r(l/2 - e/2) 
2i-^(27r)8-2T(l - e/2) ^ 




dk^ 


+ W ■ K - P(3K2 + K ■ W) 
_- 

^ (K2 + P)2((K + W)2 + P)(K2)e/2 


»1 


c^ 2 /(l -|/)x 


X 


dk. 


d2-^K + W ■ K K2 - fc2(3K2 + K ■ W) 

^(K2)^/2 [(K + 2/W)2 + W2|/(l - I/) + A:2]3 

„ 7/ ' ' I I #- / 7. I I I I / 7. - t- //. \ ! t ! I I /* ^ 

= 4(1 - iV/) 


.,7r7/2-^/2r(e/2)r(l/2-e/2) 


yJ 

dK 


X 


(K2)e/2 


2i-^(27r)8-2T(l - e/2) c 

+ W K K2 

[(K + yWy + W^y{l - y)f/^ ~ [(K + yWy + W2y(l - y)f/^ 


dy{l-y)x 

3K2 + K W 


= 4(l-iV/) — 


7r^/2 ^/^r(l/2 — e/2) 


2i-^(27r)8-2T(l - e/2) c 

+ W K K2 


dkr 


dz 


dy{l - y)d^-^K 


[(K + yzWy + y^z{l - z)W^ + W^yz{l - |/)]5/2+e/2 r(5/2) 


r(5/2 + _ ^y/ 2 - 1 _ 


3K2 + K W 


[(K + yzWy + y‘^z{l - z)W^ + W^yz{l - |/)]3/2+^/2 r(3/2) 


r(3/2 + _ 2 )^/ 2 -i 


(A6) 


Now we shift K —)■ K — i/zW. This leads to the replacements ((K ■ W)^ = K^W^/(2 — e) as far 
as integration over K is concerned) 


+ W • K ^ + -- - -+ TK^W'^y^z^ 

- 

^ + W%^z\ 

K • W ^ -yzW^. 
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Then, integrating ont K, 


{V»)seM=4(1-]V,-')- 


ir0/=-r(l/2 - £/2) -■ '' 


-u^ / dkx / dz dy 

Jo Jo 


X - 


2-^(27r)8-2^r(l - e/2)2 c 
(3|/V - yz) (r (e + I) r (1 - f)) {y^{l - z)z + (1 - y)yz)~"~^^^ 

2 r(|) 


{l-y){l-zy^^ 


+ 


3yz (r (e + I) r (2 - f)) (|/2 (^ + 2) - (^ + l)) {y\l - z)z + (1 - y)yz) " z^/'^ 

2f© 

3 (r (e - i) r (2 - f)) (?/2(i - z)z + (1 - y)yzf^‘^~" z^!'^ 

2 r(|) 

3 (r (e - I) r (3 - f)) {y^il - z)z + (1 - z^!’^ 

2 r (I) 

3y^z^{yz - 1) (T (e + |) T (l - f)) (y^(l - + (1 - y)yz)~"~^^^ z^l‘^ 

2 r(f) 


(AS) 


To leading order in e, we can take only the (1 — term in the above integrand for its z 

dependence and set z = 1 elsewhere (which prodnces 2/e for the integral over z). This agrees with 
nnmerically evalnating the y and 2 ; integrals. We thns get, 


{JyJy)sE{uj) = j dK c dy y{6y-5)^l^^ 


(A9) 


3. Vertex correction to (JyJy) 

As in Appendix A 2, here too we will freely take the limit of vanishing velocities associated with 
Eq. (3.13). The case of hnite velocities will be presented in Appendix C. 

This correction is then given by Fig. 2(c): 


^JyJy) verti,^') 2 ^Aj 


(27r)2 


d^-^K ^ 
-Tr 

(27r)2-^ 


Id-r 


J-'- /d-i' 


K2 


U2 

Ky 


(K,k,W 


(K + W)+7d_ifc, 

(K + W)2 + A:2 

(AlO) 
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We again first compute the current {Jy) vertex, given by Fig. 3(b): 
H 3 (K,k,W) 


2 „. 
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Nf — 
J j=i 




d\ d^-^Q . 


X i—t 


(27r)2 (27r)2-^ 

.^r-(K + Q) - ^d-l{qy + ky) 


*7d-iX 


(K + Q)2 + {qy + kyf 


w .N r ■ (K + Q + W) — 7rf_i(gy + /cj^) . 1 

(K + Q + W)2 + (g^ + y2 'T''i-lQ2 + c2q2 


.2,,. ^c -1 


= -2*^ V (rW^) [\x f "" 

Jo Jo 


Nf — 

I j=l 


dy 




+ 2r ■ (K + Q)7c;_i(g, 


(27r)2 (27r)2- 

, , , r-W 7 ,_i(g, + A:,)-r-(K + Q)r-w 


(K + Q) ■ (K + Q + W) + (gj, + kyf 


X 


X 


(Q + (x + 2 /)K + |/W)2 + W^yil -y) + {l-{x + |/))(K2(x + g) + 2K • Wg) + 


+ (x + y){qy + kyf + (1 - (x + y)Wc^ 


-3 


ld -1 






E 


r>l —X 


Jt-J' 


r-'r- 




rf|/- 


dqy 


cNf{27r)^-^T{l-e/2) ^ Jq Jq ^(1 - (x + i/)) 

(A:, + - (K(x + g - 1) + W(g - 1)) • (K(a; + y - 1) + Wy) + 

+ 2r ■ (K(l - (x + y)) - Wy)j,_f{ky + g,) - F • K(1 - {x + |/))r • W + F ■ W 7 ,_i(A:, + g,) 


X 


X F(3/2 + e/2)F(l - - F(l/2 + e/2)F(2 - j 7 ^. 1 , 

where we again integrated out q^ and then sent c —)■ 0 in the last step of the above, and 
Ai = W2|/(l - ^) + (1 - (a; + y)){^\x + |/) + 2W ■ Ky) + {x + y){ky + qy)^ 
Proceeding, 

A/2_1 

1 


(All) 


(K,k,W) = -z- 




iV 2 -l 

E 


'j-J ^J'' 


'•I pl — X 

dx I dy- 


cA/(27r)^ T(1 - e/2) ^ Jq Jq ^(1 - (x + y)){x + y) 

- (K(x + g - 1) + W(g - 1)) ■ (K(x + g - 1) + Wg) - F ■ K(1 - (x + g))F • W 


X 


X r(i-./2)r(i + ./2) _ _ r(i-^./2) ^ 


Ad+^/ 2 ) ^72 2(x + g) 

where 

A 2 = W2g(l - 2 /) + (1 - (a; + g))(K2(x + g) + 2W • Kg). 
An important feature of the above computation is that because 

2(x + g) - 1 


7/2 


'‘I pl—X 

dx I dy 


'0 


Cl (x + yf/^^1 - (x + y) 


= 0 , 


(A12) 

(A13) 

(A14) 


30 




























the coefficient of the 1/e pole (i.e. the coefficient of r(e/2)) in H vanishes when e —)■ 0. This 
eventnally leads to the lack of a 1/e pole in {JyJy)yeTt, and hence the correction to scaling of (JyJy) 
arises solely from the self-energy graphs. 

Taking the expression for the cnrrent vertex and inserting it into the one for {JyJy)vert, we get 


^JyJy) vert ) 


dk^. d^‘^'K dz dx 


'*1 — X 


dy 


c(27r)8-2^r(l - e/2) 

- (K(a: + y-l)+ W{y - 1)) ■ (K(a: + y - 1) + Wy) 


0 Jo y/(l - (x-F|/))((x + |/) 

/ 1 K ■ (K + W) 




1/2 


.3/2 

^3 


+ 


K^W^l - (x + y)) 
aT 2 


r(i- 6 / 2 ) 

2{x + y) 


r(l-e/2 )r(l + e/2) r(e/2) 


(l-l-e/ 2 ) 




r( 2 -£/ 2 )- 


a 1/2 
^3 


K ■ (K + W) 

^3 


(A15) 


where now 

A 2 = W2i/(1 - ^) + (1 - (x + y)){K\x + y) + 2W- Ky), 
As = (K + zWf + z{l - z)W^. 

We combine denominators nsing 


T{s + b) 


da-, 




b-l 


A|A^ r(s)r(6) io [aAa + (1 - 

and the denominator sqnare completion is 

'W{ay{l - (x + y)) + (1 - a)z) 


(A16) 


(A17) 


aAa + (1 - a)A 3 = {a{x -\-y){l-{x + y)) + (1 - a)) 


a(x -f i/)(l - (x + y)) + (1 - a) 




Dehning 


(ay(l - (x + y)) + (1 - a)zY 

a(x + y)(l- (x + y)) -1- (1 - a) 


+ a(l - y)y+ (1 - a)z 


KJ, 

(A18) 


/i = a(x -hy)(l- (x + y)) + (1 - a), 

/ = - (^ + y)) + (1 - 

/2 = a(l - y)y + (1 - a)z - /l/^ 

we can process the numerators and write down the final expression 


(J^pJ^p)vert(^) 


- 1 ) 1 - 


f*l — X 


c(27r)8-2W(l - e/2)‘ 


u 


2 e 


dkr 


da dz dx 


(A19) 

dy{ + ^2 -l- Ts -|- , 

(A20) 
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where 


--S 

n o p e ^ 


T = yi - - 

^ ^J-{x-Yy- l){x + y)TT V 2J\e-2 


fJ,T{2--jT{e+-] {2f\e-A){x + y-iy 


- /(e - 4:){x + 3y - 2){x + y - 1) - 2xy + X + {y - l){y{e - 4) + 1) - (1 - (x + y)){e - 1) ) + 

+ (/ - i)//fr (i - 0 r (f +1) (/(^ + 9 -1) - v){f{x + 9 -1) - 9 + 1)+ 


/^(a: + !/-l)"r 3-- r 6-- 


r, = 


a-'VrV,"" 


=r 


vT^ 

X {fix + y-l)-y + l)- f 2 {x + y- 1)^1 ^ - 0 j , 


(i-|) (-/ir'(i-|)r(€ + 0 (/(x + !,-i)-9)x 


T, = 


\/l - 002 1/^2 2 ^^ . 2 


e\ ‘ 


-= VW-x-9+11x + 9)^/- "^^-2J r^.-||W.-2) + 9(.-2) + l)|/^x 
X (/i — 2/ie) + //i(2e — 1) + / 2 (e — 2) ), 


r4 = 


/2 - ^r(l-) (x(e-2)+|/(e-2) + l). (A21) 

4\/l — a-^/ 7r(—X — 1 / + l)(a; + y)^/'^ ^ 2/ \ 2 


This multidimensional integral over four parameters is finite in the limit of e —?■ 0 and can be done 
numerically. We hrst integrate over x and y: The resulting function of a and +; has integrable 
singularities in the limits of a —?■ 1 and a —?■ 0 which can be handled by numerical integration using 
an adaptive grid. The final result is 


{JyJy)verti^') I d/c. 


- 1 ), , 1 - 2 . 


327r4c 

where cio ~ 1-1 is a hnite numerical constant. 


-u 


oo + I) 


(A22) 


Appendix B: Free Energy Computations 

As in the previous appendix, we will freely take the limit of vanishing velocities associated with 
Eq. (3.13) here as well to compute the correction to the fermion free energy. 


Ffb = -Tr 


'n^-1 

f (Pqd^ 1 

2_^T T 

i=i 

J (27r)^“^ Q2 + c2 q 

LJq 




[(n(g,T)-n(g,0)) + n(g,0)]. (Bl) 
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Where 11 (g, T) is the fermion RPA bubble at external momentum and frequency given by q evalu¬ 
ated at temperature T. As described in the main text, we evaluate the hnite temperature part of 
the bubble at n = 0 to renormalize the fermion free energy and the zero temperature part at n 7 ^ 0 
to renormalize the boson free energy. To evaluate the frequency summations, we use the following 
zeta-function regularization identities: 






1 



27 —^ 

{27ry 


■'1 — 5 


( 2 ir) 


:C(o), 




(B2) 


Where ojq is a bosonic Matsubara frequency and ook is a fermionic Matsubara frequency. We then 
have 


n(g,T)-n(g,o) = -V 
= / dk^ 


dkx f dkyd} 


dkyd} 

3-e 


( 2 vr: 


2 vr , 


(27r)2-^ 
duk 






duk 

2ti 


Tr [i^d-iG{k)i^d-iG{k -h q)] 


‘^k 


271 


Q^ + ql+co; 


(K^ + ky + ci;^)((K -|- Q )2 -I- [ky -\- qy^ -|- [uk + 
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1 


ky uj'l (K -1- Q)^ -|- [ky qy^ -\- [uk -\- ojq)^ 


(B3) 


The last two terms in the square brackets yield identical contributions, because the q in the last 
term can be shifted out. Thus, 


n(g,T)-n(g,0) 
= / dk^ 


dkyd^ ""K I ^ y-v 


( 2 vr: 


3-e 




-I- / dk^ 


dkyd^ I ^ y-v 


( 2 vr: 


3-e 


‘^k 


r (k^ 

I 271 

doJk _ 

27r ] K^ + kl + ul' 


^^ + ql+ujl 


(K2 + kl + a;2)((K + Q)2 + [ky + g,)2 + (a;, + a;,)^) 
1 


+ 


(B4) 


We evaluate the second term to leading order in e in the above using dimensional regularization 
for the momentum integral and zeta function regularization for the frequency sum: 


^9^ 



dkyd^ 

(2vr)3-^ I ^ 

\ <^k 


duk\ 1 

^ ) K 2 + + OjI 
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dkyd^ V 1 
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(B5) 
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where we have used the fact that scaleless integrals vanish in dimensional regularization. Thus, 


n(g,r)-n(g,o) 



dkyd^ 

(27r)^“^ 




8g^T^-^ In 2 
(27r)2 




_+ qj + _ 

(K2 + ky + a;^)((K + Q)2 + (/c^ + gy)2 -|- + Ug)'^) 

(B 6 ) 


To evaluate the hrst term, we introduce a Feynman parameter y to combine the denominators. 
Doing the k momentum integral and Uk frequency summation (integral for the T = 0 part), we 
have, to leading order in e 


n(g,T)-n(g,0) = 


Stt^ 


dk^ 




e,g‘T'-—\n2 


(27r) 


(B7) 


Where q = (Q^ + We determine the following asymptotic expansion numerically 


dy t{y,f,ujg,e) = 27r - T 
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(B8) 


Simple power counting dictates that the higher terms in the above asymptotic expansion can’t 
produce any log UV divergences in the hnal two loop graph because they fall off too fast in q. 
Thus, retaining only terms that will survive and contribute to the pole in the hnal two-loop integral, 

n(,, T) - n(,, 0) = -|1 y di, (^48C(3)r»->j , ( 39 ) 

We evaluate n(g, 0) using dimensional regularization at hnite v to get, to leading order in e: 
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r 
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V Svre 
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(BIO) 


where the last integral was performed using Feynman parameterization. Inserting the expressions 
for n(g,T) — n(g, 0) and n(g, 0) into Eq. (Bl), we get, using dimensional regularization for the q 
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momentum integrals, 
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Where we integrated out and then sent c —)■ 0 in the last step of the above. Doing the remaining 
integrals first over qy and then Q, we get, to leading order in e 


( 1 ) - l)T 
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The other part gives 

- 1i _ n(g,0) _ g'^ r 
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(B13) 


Integrating first over Q and then over q using the dimensional regularization, we get, to leading 
order in e 
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( 2 ) 

fb 
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|-(3-2e) 




(Nc - 1 ) 


ShOuc^e 
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Appendix C: Finite v and c 


In this appendix, we describe the breakdown of the results derived in the previous appendices 
when we do not have u, c —?■ 0. We illustrate this by first computing the self energy correction to 
(JyJy) for finite v and c; similar problems occur in the computations of the fermion free energy. 
The fermion self energy for u, c 7 ^ 0 is given by [36] 


Si(K,k) = - 
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2-e 


I^T(e/2) 




N'^-l 
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/ Jo 


V- r K 
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I c2£;2(k) 

c^+x{l+v‘^-c^) 


^/2 (c2 + x(l+u2_c2))l/2’ 

(Cl) 


We can ignore the term with the prefactor of F in the numerator of the integrand in Eq. (Cl); 
since u 7 ^ 0 £i(k) and e 3 (k) can be taken to be independent variables of integration over k space 
via the coordinate transformation —)■ dside^/ (2u). This term then only produces contributions 
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to {JyJy)sE that are odd in £3 and hence vanish under integration over £ 3 . Thus dropping this 
term, we have 




16(1 - iV2)7r2-^/2r(e/2)c/V 
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deide^i 2 -^ 
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d^-^K- 
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(C2) 


evaluating this as in Appendix A 2 gives the singular contribution 


{Jyjy) Se{^) ~ 



,2^2 


C 8 


u'^{c^ + x{l + n 2 


_ \ ( 1 -x)V 2 ( 1 +^ 2 ) / 1 X 

— c2))^/2 y647r^e/ ’ 
(C3) 


where the crossover function Ki{x, e) ~ (1 + for x <C 1 . 

The singular part of the 1-loop current vertex at hnite v and c is most easily derived from the 
Ward identity; 
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(C4) 


Inserting this into Eq. (3.16) gives, for the singular part of the two-loop vertex correction to (JyJy), 
via a computation very similar to that for the self energy correction (an additional prefactor of 2 
has to be inserted to account for both the poles associated with vertex corrections to each of the 
two current vertices in the graph). 
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where again the crossover function K 2 {x, e) ~ (1 -|- x) ^^2 for x 1 . 
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Appendix D: Boltzmann Equation Computations 


1. Collisionless conductivity in d = 3 — e 


We can diagonalize the Hamiltonian corresponding to the free fermion part of Eq. (2.3) as 


4 Nf 




= E E E / n H ■ K + K) 

n=l a=l j = l J 

Nf - - - 

= E E E E / K)?„,„(k, K)A„,,j,„(k, K) 


n=l <7=1 j=l m=± ' 


(Dl) 


with the particle-hole symmetric dispersions K) = m (K^ -|- The physical current 

density becomes 

J = EEEE / 9 ?? +n. (d 2 ) 


n=l cr=l j=l m=di ' 


VK2-F£„(k)^ 


where £n(k) = ■ k and J 2 contains particle-hole terms A!)_A_, A1A+ that are unimportant for 

transport in the low frequency regime of cn T [54, 55]. Defining the distribution functions 


/n,m(k,K,t) = (A^j-^(k,K,t)A„,^j,m(k,K,f)), 

we have the collisionless kinetic equation in the presence of an applied electric field 


(D3) 




^ + ™E' ^ I /n,m(k,K.t) = 0, 


(D4) 


with the frequency-domain solution to linear order in E 

m£„(k) 


/n,m(k, K, cj) 27r(5(cn)?T.j'(^^ ni(k, K))-|-v^'E(cj) - 


1/T 


v/K2-t-£„(k)2 -iu: + 0+ 
Inserting this into the expression for J, we obtain the collisionless conductivity 


'>/K»,m(k, K))(l-n/(^„,„(k, K))). 
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= 2N,NjVTTV^ j ~ = Re|a„((.,)]. (D6) 
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2. Derviation of the fermion collision integral 


We derive the following expressions for the different components of the fermion self energies in 
the Keldysh formalism. For a fermion at hot spot given by (£, +), we get for the hrst diagram for 
Sj in Fig. 5 

= ig‘^T^pTp^,D^{x,x')G^-{x,x') = ?,i5„„ig‘^D^{x,x')G^-{x,x’). (D7) 

We use that for products of Wigner transforms 


D^(x,x')G^ {x,x') ^'^Dl^{x,p-q)GS‘ {x,q), 


(D8) 


and plug in the representation of the Keldysh propagator in terms of the distribution function to 
get 


Y:^^J;^\x,p) =36aa'g'^^Fb{x,p-q)i [D^{x,p-q) - D^{x,p- q)] G^ {x,q). (D9) 

q 

For the collision integral on the right hand side of (5.16) we need twice the imaginary part of this 
expression. Using 


2Im[G'f {x, g)] = - [G^ {x, q) - Gf (x, g)] , 
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(DIO) 


we get, 


2Im 
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= 25^^/Im 


(Dll) 
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where we have used spatial translational invariance and also have kept the external fermion on 
shell. Likewise, for the second diagram for in Fig. 5 contributing to we have 


2Im|E“:P'(i,p)] 

= - G"~(x,q)]lDllix,p-q) - D^{x,p-q)] 


= -3(5a<7'5'^ 
= 25^^/Im 


f (fc[ 1 

J 2,71 4cnp_q 

E“+^^^(f,p,e+(p)) 


((5(e+(p) - (q) -Wp.q) - (5(e+(p) - (q) +a;p_q)) Fj (t, q) 


(D12) 


For the diagrams in Fig. 5 contributing to the hrst gives 


^ Ff{x, q)Fb{x, p-q) [G^” (x, q) - Gq ^"(x, g)] [D^{x, p- q) - D^{x,p- q)] 


r X 
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(D13) 


The second and third combined yield 


= -^5„^ig‘^^[G^-{x,q)D^{x,p-q) + Gf-D^{x,p-q)] 
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+ c.c. 
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Combining the above expressions gives the collision integral for fermions of any spin at the hot 
spot given by (£, +) 


= 3g^ f 


4^p—q 


^5(e/(p)-e^ (q)-a;p_q) 1 + FJ (t, q)F 5 (t, p - q, a;p_q) - (t, q)Fj+(t, p)- 


Ff^{'t,p)Fb{t,p-q,ujp-q) - (5(e|(p) - (q)+u;p_q) 1 + Fj (t, q)Fb(f, p - q,-Wp.q) 


Ff {t, (l)Ff^{t, p) - p)Fb{t, p - q, -Wp-q) 


?e+ 


(D15) 


The collision integral for the hot spot given by {i, -) is given by simply interchanging + f-)- — in 
the above. With the alternate parameterization (5.15) of the distribution functions, and relabeling 
of momenta q -H- p — q we may rewrite this for any hot spot as 

= 3c/^ j ^<5(e^(p) - e^(p - q) - Wq) p)(l - p - q)) - p - q)/b(t, q, u;q) + 

+ /^^(Cp)/fe(Cq,Wq) - <5(e^(p) - ef{p - q) + Wq) /^^(f,p)(l - ff{t,p~ q))- 

-/^^(Cp-q)/6(Cq,-Wq) + /J^(Cp)/b(^,q,-Wq) )• (Die) 


3. Derivation of the boson collision integral 


We begin with the retarded component of the boson self energy in the Keldysh formalism. The 
sum of the two diagrams for this component in Fig. 5 gives 

Sf (x, x') = -ig‘^ ^ x')Gf~{x\ x) + G^~{x, x')Gq^^{x', x) + {+^ -)] . (D17) 

i 


Wigner transforming this gives 

2Im[Sf (x, q)] = -g q) (Gf+(x, k + q)- Gf+(x, k + q)) {Gf {x, k) - G^ {x, k)) 

l k 

+ Ff{x, k) (Gq^^{x, k) - Go^+(x, k)) [Gq^~{x, k + q) - Gq^~{x, fc + g)) + (+ -H- 


= -9 




27r L 


S (e^ (k) +a;q-e/(k + q)) k + q)- 


e •' 

— 6 {e\ (k) + (Uq — (k + q)) Fj~^{t, k) + (+ -H- —) 

= 2Im[Sf(t,q,a;q)], (D18) 
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where we have used spatial translational invariance and also have kept the external boson on shell. 
For the Keldysh component of the boson self energy, the second diagram in Fig. 5 gives 


= ^ k + q)- k + q)) Ff{x, k) {G^~{x, k) - Gf~(x, k)) + (+ O 


l k 


i 


(f\i 


2ti 


[(5 (cjq + (k) — e^(k + q)) k + q)^^^ (t, k) + (+ -H- —)] 


(D19) 


The hrst and third diagrams for the Keldysh component when combined give 




K{l+3) 


[[Gq^~^{x, k + q)GQ^^{x, fc)) + {Gq^^{x, k + q)G^ {x, fc)) + (+ -H- —)] 
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d^k f du 


dTT^ / 271 


+ C.C. + (+ O -) 


(Uq + a; — (k + q) + i0+ OJ — (k) — i0+ 

= g'^ ^ f [d (wq + (k) - e/(k + q)) + (+ f-)- - 

Thus we obtain the boson collision integral: 

Fk] {t, q) = ^ <5 (e7 (k) + Wq - e+(k + q)) 

^ - Ff~{t,'k)Fh{t,q,u^) + Ff-^{t,k + q)Ff-{t,k) + Fj^{t,k +q)Fb{t,q,uj^)) + (+ O -) 


(D20) 


(D21) 


or 


d^kr 


fb]{t, q) = V ^ (e^ (k) + Wq - e+(k + q)) k + q)(l - /^ (t, k)) + 


+ k + q)fbit, q, a;q) - fj (t, k)/b(t, q, Wq)) + (+ ^ -) 


(D22) 


4. Solution of the boson kinetic equation 

Inserting the parametrization Eq. (5.22) for the fermion / functions into the boson collision 
integral and also parameterizing ff, in the frequency domain as 

/fe(a;,q,a;q) = 27rd(a;)nfe(a;q) + n(a;, q, Wq), (D23) 
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where u is linear in E. Changing variables to (k) = • k and P 2 — = v/ ■ k, the 

boson collision integral (5.20) becomes 



Always integrating out P 2 expression, and keeping only terms up to linear order in E, we 

get, in the frequency domain, using the boson kinetic equation Eq. (5.17) 


2(-za; + 0+)M(a;,q,u;q) = 



+ E(a;) • Wq, ■ q)nb(a;q) + E(w) ■ ■ q) - E(a;) ■ d\p{,Uq, ■ q)- 


c{p{,ujq,^i • q)(2vr(5(a;)) + ai(pi,a;q, • q)(27r(5(a;)) , 


(D25) 


where 


a(p^i, Wq, ■ q) = (ufipi + Wq) - nf{p{)) + (nf(p{ + • q) - nf{pi + • q - a;q)). 



b^(p^i,a;q, v7 ■ q) = v+n/(p^i + a;q)(l - nf{p{ + Uq)Mp{ + Wq) - vjnf{p{){l - nf{p{))p{pi) + 
+ v7 nf{p{ + v 7 • q)(l - nf{p{ + • q))cp{p{ + ■ q)- 

- v+n/(pi + v7 • q - u;q)(l - nf{pi + ■ q - a;q))(^(pi + • q - Wq), 

bi(pi,a;q, v7 ■ q) = v+n/(p^i + a;q)(l - nf{p{ + uj^))p{p{ + ujci) + 

+ nf{p{ + vj ■ q)(l - nj{p{ + ■ q))^{p[ + ■ q), 

d^(p^i,a;q,v7 ■ q) = wjnf{p{+ uj^)nf{p{){l - nf{p{))p{p{) + 

+ + ‘^q)(l - nf{p[ + ^q))p{p[ + C^q) 

+ • q)nf{p\ + v7 ■ q - a;q)(l - nf{p[ + • q - uj^))p{p[ + ■ q - a;q) + 

+ nf{p[ + v7 ■ q - a;q)n/(p^ + ■ q)(l - nf{p[ + ■ q))^{p[ + ■ q), 

c(/i, Wq, v7 ■ q) = nf{j)[ + UJq)nf{j)[) + nf{j)[ + ■ q)n/(p^i + ■ q - Wq), 

ai(pi, Wq, v7 ■ q) = nf{p[ + Wq) + nf{p[ + • q), 

/ dp[ [ai(pi,a;q, • q) - c(pi,a;q, ■ q)] = 2 / dp\nf{p[ + a;q)(l - nf{p[)) = 2u^nb{u^). 


(D26) 


42 






Since each term in the b^, b^, terms results in a convergent integral over the ■ q s can 
be shifted out. Then, since = 0, the contribution from the b^, b^ and terms vanishes. 

Then, 


2{—iuj + 0^)m(ci;, q, cuq) = —8uci—u{u, q, Wq). (D27) 

TIV 

Since this has to hold for all values of uj, we can only have u = 0. Hence, the boson collision 
integral is trivially solved by the thermal Bose distribution and the bosons do not respond to the 
applied electric held in our approximation. It is also easily seen using the identity 


nf{x){l - nf{x - y)) + nb{y){nf{x) - nf{x - y)) = 0, (D28) 

that the thermal Fermi distribution Uf nullihes the fermion collision integral in the absence of an 
applied electric held if the thermal Bose distribution nj, is used for the bosons, as it should. 
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